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Two-dimensional (2D) particulate aggregates formed due to competing interactions exhibit a range 
of non-equilibrium steady state morphologies from hnite-size compact crystalline structures to non¬ 
compact string-like conformations. We report a transition in heterogeneous microscopic dynamics 
across this morphological hierarchy as a function of decreasing long-range repulsion relative to 
short-range attraction at a constant low density and temperature. Following a very slow cooling 
protocol to form steady state aggregates, we show that geometric frustration inherent to competing 
interactions assures non-ergodicity of the system, which in turn results in long-time sub diffusive 
relaxation of the same. Analysing individual particle trajectories generated by molecular dynamics, 
we identify caging dynamics of particles in compact clusters in contrast to the bonding scenario for 
non-compact ones. Finally, by monitoring temperature dependence, we present a generic relation 
between diffusivity and structural randomness of the aggregates, irrespective of their thermodynamic 
equilibrium. 


I. INTRODUCTION 

Relating the morphology and dynamics of non¬ 
equilibrium systems has remained a challenge despite 
concerted theoretical and experimental efforts over 
the years. Systems driven out of equilibrium, such as 
glassy and granular materials, exhibit dynamic structural 
order mi] which often extends over different length and 
time scales without any clear separation between them. 
The effect of geometry on dynamics mil] is also expected 
in systems confined by external fields or porous media. 
Transport within crowded environments such as living 
cells [9] presents yet another example where spatial struc¬ 
ture and dynamics are surrogate to each other. Here, we 
present a detailed study of both global and local dynam¬ 
ics of a generic non-equilibrium pattern forming particu¬ 
late system, namely, 2D aggregates. Aggregates of parti¬ 
cles formed due to competing interactions enjoy a broad 
range of structural randomness (or order) [TOHE] inter¬ 
mediate between liquids and crystalline solids. By fol¬ 
lowing the steady state trajectories simulated by molec¬ 
ular dynamics, we establish a phenomenological corre¬ 
spondence between the local morphology and dynamics 
of 2D aggregates. 

Experimentally, the aggregates are usually formed by 
quenching a high-temperature equilibrium liquid [laiii] 
to a much lower temperature. While such a dynamical 
preparation protocol is useful in practice, identification of 
the governing dynamics and a theoretical understanding 
of the processes becomes difficult. Different mechanisms 
such as diffusion-limited [20] and reaction-limited aggre¬ 
gation m have been reported. To describe the phase be¬ 
haviour of such systems, both equilibrium routes [2^1^ 
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and kinetic viewpoint [2H426] have been thoroughly ex¬ 
plored. Stemming from a fast quench, the glass transition 
scenario has also been conjectured and quite thoroughly 
investigated EilllHl. 

In the present work, a very slow cooling protocol is 
adopted for aggregate formation and the final temper¬ 
ature is maintained over a long observation time upon 
reaching the desired value. By monitoring the evolu¬ 
tion of energy fluctuations during cooling, we show that 
frustration innate to competing interactions is the key 
factor for the non-equilibrium dynamics of aggregates. 
The steady state dynamics is non-ergodic, sub-diffusive 
and heterogeneous for all aggregate morphologies. While 
phase separation affects this overall dynamics, local slow 
dynamics results from the reduced thermal fluctuations 
at final ‘low’ temperatures. We find that individual par¬ 
ticle trajectories and their spatio-temporal correlation re¬ 
veal a continuous transition in microscopic dynamics in 
concert with the changing shape of aggregates. Finally, 
we demonstrate the scaling of the average dynamics spec¬ 
ified by scaled diffusivity of aggregates with respect to 
their local geometry quantified by the excess entropy. 
This relation points to a fundamental connection between 
average dynamics and structural randomness of a system 
irrespective of its thermal equilibrium. 

The organisation of the paper is as follows: we intro¬ 
duce the specific form of competing interactions used in 
the study in Sec.II and briefly describe the variations of 
morphologies obtained by tuning the competition. The 
details of cooling protocol and molecular dynamics simu¬ 
lation are charted out in Sec.III. In Sec.IV, we elaborate 
the main results on dynamics in three parts, (A) non- 
ergodicity, (B) global dynamics and (C) particle dynam¬ 
ics. A discussion and further implications of these results 
are presented in Sec.V followed by a brief summary. 



2 


II. VARIATIONS IN MORPHOLOGY FROM 
COMPETING INTERACTIONS 


Competing interactions are necessary for aggregate for¬ 
mation Eiisn] and can occur in diverse physical settings 
[ 3 TII 34 ] . each having its own set of control parameters 
[ 35 ] . Our choice of model system is a reliable represen¬ 
tative of globular proteins [36] and colloidal systems m\ 
with suitable polymer coating. In this model, the short- 
range attraction and long-range repulsion are realised by 
a generalised Lennard-Jones potential [38] and a pair¬ 
wise potential of Yukawa form, respectively, 

(/.5A=4e[(cT/r)2«-(cT/r)“] (1) 

(t>LR = (^cr/r) exp(r/C (2) 

Setting length and energy scales by a and e respectively, 
the range of attraction is controlled by a. The attrac¬ 
tion range is 0.2cr for our choice a = IS. The same ex¬ 
tends to 2.5cr for the usual Lennard-Jones potential with 
a = 6. We note that the thermodynamic properties of 
a particulate system with bare attraction is independent 
of a ioT a > 18 [38]. The Yukawa part implicitly mod¬ 
els the long-range repulsive effect of surrounding media. 
The strength of repulsion A and the screening length ^ 
are measured in e and a units respectively. The effec¬ 
tive centro-symmetric interaction, (j) = (t)sA^(t>LRi has a 
positive energy barrier at finite distance separated from 
the steep hard core repulsion by an attractive minimum. 
The height and extent of the barrier can be controlled 
by A and In a previous study [39], we have shown 
that this system exhibits a structural hierarchy at a con¬ 
stant density and temperature for different realisations 
of competition controlled by a, A and Aggregate mor¬ 
phology shows a continuous transition from finite-size 
non-compact to compact to percolated gel structures as 
a function of decreasing repulsion against attraction. 

In this study, we focus on a restricted parameter space 
of only finite-size clusters and investigate the dynamics 
across compact to non-compact morphologies. Specif¬ 
ically, we present the detailed analysis of microscopic 
dynamics for two representative cases, ^ = 0.5 and 0.8, 
showing compact and non-compact aggregates respec¬ 
tively for fixed A = 4.0. We mention that the strength 
of repulsion equals that of attraction for this choice of 
A. The form of the effective pair-wise interaction (j) is 
shown in Figj^a). Note the dominant attraction well 
for ^ = 0.5 which is complemented by a small repulsive 
barrier slowly decaying at large r. The resulting con¬ 
figuration consists of perfect crystalline islands of expo¬ 
nentially distributed size. Percolated structures appear 
for lower values of ^ which we keep out of the present 
study. When ^ is raised to 0.8 keeping A constant, the 
repulsive barrier increases significantly and the attractive 
minimum close to the particle core (r = 0) becomes com¬ 
parable with the repulsive minimum far from the core. 
Finite-size one-dimensional arrangements of particles are 
observed under such strong influence of repulsion. Upon 
further increase of the size of non-compact chain-like 
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FIG. 1. (color online) (a) The form of effective interaction 
0(r) is drawn for two different values of ^ = 0.5 and 0.8 seting 
A = 4.0 and o = 18. Please note that the global attractive 
minimum becomes local with increasing Visual represen¬ 
tation of particle dynamics for these two ^ values are shown 
in (b) compact (^ = 0.5) and (c) non-compact (^ = 0.8) struc¬ 
tures, respectively, at T = 0.05. Only a part of the simulation 
box is shown for clarity. The scale used for both (b) V (c) 
is same and shown in (c). Each red filled circle denotes an 
instantaneous particle position. Blue background consists of 
dots representing the positions of the same particles for 100 
earlier instants stored in a regular interval of 500 simulation 
steps, (e) V (f) are trajectories of a single particle belong¬ 
ing to (c) & (d) respectively, colour coded by time. Using 
the same scale (scale bar = 0.5cr), it becomes evident that the 
particle trajectory (d) in compact clusters is more restricted 
and jagged than that in (e) non-compact clusters. 


aggregates becomes smaller and the overall randomness 
in the configuration increases further. 

A part of the simulation box is presented in 
FigQb)&(c) showing both of the limiting morphologies, 
compact and non-compact, respectively. Red Ailed circles 
in these figures represent instantaneous positions of par¬ 
ticles. Positions of the same particles for 100 previous 
instants are plotted in the background with blue dots. 
All the collated configurations are separated from each 
other by 500 molecular dynamics steps. The slow dy¬ 
namics experienced by the particles is already evident 
from these visual representations. Another important 
observation is the consistency of collective dynamics of 
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particles. Particles forming a cluster undergo collective 
translation and/or rotation; even bending to change con¬ 
formations. However, fragmentation and reformation of 
aggregates are seldom observed. In other words, the sys¬ 
tem has reached its steady state at least in terms of the 
size of clusters. Trajectories of a single particle chosen 
arbitrarily from the compact and non-compact clusters 
are drawn in Figj^d)&(e) respectively with colours ac¬ 
cording to their time evolution. The spatial extent of the 
trajectory is restricted for compact aggregates compared 
to that of a non-compact one. The former appears to be 
more heterogeneous than the latter from visual inspec¬ 
tion. We analyse these trajectories in terms of several 
measures of correlation to understand the microscopic 
dynamics of the systems. Before discussing those results, 
we briefly describe the cooling protocol observed in this 
study and comment on the energy fluctuation during and 
after cooling in the next section. 


III. COOLING PROTOCOL 

We consider a 2D system of interacting mono-disperse 
particles at fixed density p = 0.4. The initial configuration 
is a random non-overlapping arrangement oi N = 56000 
particles of unit mass in a 376.0 <j x 372.0cr box with pe¬ 
riodic boundary conditions along all directions. Parti¬ 
cle trajectories in canonical ensemble are generated by 
following molecular dynamics (MD) at constant number- 
area-temperature {NAT) as implemented in LAMMPS 
simulation package m- Temperature is measured in e 
units with the Boltzmann constant, ks = ^ and main¬ 
tained by a Langevin thermostat m- While the parti¬ 
cles experience forces from others due to pair-wise inter¬ 
action 0, they are also subject to damping from implicit 
surroundings and a random force of Gaussian variance. 
The motion of i-th particle with position vector can 
then be expressed as follows, 

ri = -^V4>{r)-uri + Ci (3) 

considering the force due to interaction (/)(r) and fric¬ 
tional drag uTi from implicit media experienced by the 
particle. Q is a random force with zero mean and Gaus¬ 
sian variance, {(i{to)(j{tto)) = 2kB{T/u)SijS{t). Nu¬ 
merical integration is performed by using velocity Verlet 
algorithm with time steps 6t = 10~^r where r = aj 
is the unit of time for unit mass. The system is first 
equilibrated at temperature = 1.0 to form a simple liq¬ 
uid. We mention that the system goes from ballistic to 
diffusive regime in time slightly longer than r and this 
feature holds over a range of temperature up to T ^0.2. 
The temperature of the equilibrated liquid is then lin¬ 
early ramped down over 9.5x10^ MD steps to final tem¬ 
perature Tf = 0.05 (Figj^a)). The temperature of the 
system thus drops by only 10“^e over unit time r. 

Instantaneous temperature or kinetic energy per par¬ 
ticle Ck is shown in Figj^a) during linear cooling. The 
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FIG. 2. (color online) (a) Linear cooling denoted by decreas¬ 
ing average instantaneous kinetic energy tk {solid blue line) 
and corresponding variation of average potential energy Cp 
{solid purple line) are plotted. Please note the sharp decrease 
in Cp near Ck = 0.2 which is close to the gas-liquid critical 
temperature for a system with bare attractive potential, (b) 
Non-exponential decay of normalised Cp {see text) at fixed final 
temperature T = 0.05 is shown for different values of ^ = 0.5 
(red), 0.6 {light blue), 0.7 {dark blue) and 0.8 {green). 


potential energy per particle e^, following e/c, shows an 
initial linear decay but drops abruptly around T ^0.2. 
We note that the critical temperature of the bare attrac¬ 
tive potential is Tc = 0.18 ± 0.01 [42]. This sharp drop 
of Cp close to Tc is indicative of the effect of phase sep¬ 
aration process on the system. After reaching a desired 
final temperature, it is maintained by the same thermo¬ 
stat over a period of lO^r. This is our observation pe¬ 
riod when the trajectories are saved for further analysis. 
Most importantly, we note that the Cp continues to decay 
and does not reach a steady state value over these three 
decades of time. This decay can be fitted well with a 
stretched exponential or a double-exponential function. 
In Fig§b), we show this non-exponential decay for dif¬ 
ferent values of We have plotted the scaled potential 
energy, Ip = ep{t)/ep{A) to compare its evolution for dif¬ 
ferent Without getting into the details of the fitting 
parameters, we would like to emphasise that our model 
system, although cooled very slowly, remains out of equi¬ 
librium over a considerable period of time. 


IV. RESULTS 
A. Non-ergodicity 

Unusual time evolution of potential energy at constant 
T prompts us to study its fluctuations further. We com¬ 
pute the energy fluetuation metrie [43[ |44| defined as fol- 
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lows: 


^lE{t) = l/Nj2{ei{t)-{e{t))f (4) 

i 

for the ensemble of simulated particle trajectories with 
average energy E over time t. F or the i-th particle, its 
average energy upto time t, ei{t) and ensemble averaged 
energy (e(t)) read as, 

ei{t) = 1/t f dt'Ei{t') (5) 

Jo 

N 

{e{t)) = l/Nj2ei{t) (6) 

i=l 

According to the ergodic hypothesis, 

ei{t^oo) = {e{t)) = {E) (7) 

Hence, an equilibrium system should fulfil the ergodic 
convergence criterion^ ^0. Following this, 

at large t, the energy fluctuation metric should behave 
as [44] 


fiB(O) Det ^ ’ 

for a system approaching equilibrium. The energy diffu¬ 
sion constant De is the rate at which the particles access 
the available conformational space. However, our obser¬ 
vation clearly deviates from this (FigJ^. The evolution 
of f^(t)/f^(0) is much slower at long time. Temperature 



remains steady with nominal fluctuations (FigJ^mset)) 
over this whole period of observation.With this firm ev¬ 
idence of non-ergodicity of the model system prepared 
by slow cooling and uninfluenced by any other external 
field, we attribute this feature entirely to the inherent 
geometric frustration of competing interactions. 


B. Average dynamics: mean square displacement 
and self-overlap statistics 



FIG. 4. (color online) (a) Mean square displacement (MSD) 
of particles is plotted against time for compact (^ = 0.5) and 
non-compact = 0.8) aggregates with fixed A = 0.4. At 
low temperature, T = 0.05, MSD shows three distinct regime, 
short-time ballistic, intermediate slow relaxation and long 
time sub-diffusive behaviour ((r^(t)) ^ , n < 1). Interme¬ 

diate slow relaxation is more prominent for compact clusters 
than the non-compact ones. Both systems behave as a simple 
liquid at high temperature T = 0.25 showing only short-time 
ballistic and long-time diffusive behaviour, (b) At low tem¬ 
perature (T = 0.05), the self-overlap function g(t) (see text) 
follows a power-law relation with time after an initial resi¬ 
dence time Tr which is larger for compact structures compared 
to the same for non-compact ones. Upon increasing tempera¬ 
ture, it behaves like a normal equilibrium liquid (at T = 0.25 
and beyond). 

We now focus on analysing the dynamics of constituent 
particles in this non-equilibrium setting. We start with 
the mean square displacement of particles over the obser¬ 
vation time. Mean square displacement (MSD) defined 
as. 


N 




(9) 


FIG. 3. (color online) The temperature fluctuation AT 
after cooling is shown in the inset. For a system approach¬ 
ing equilibrium, the normalised energy metric decays linearly 
with inverse of time (see eq^. However, the long time decay 
is evidently much slower and thus the non-ergodicity of the 
aggregate systems over the observation time is established. 


measures the spatial extent of a particle’s displacement 
from its origin after which the particle diffuses randomly 
due to thermal fluctuations. Time dependence of MSD 
is generically expressed as (r‘^(t)) ^ E [45] where free 
diffusion is observed for n = 1 and n < 1 indicates sub¬ 
diffusion. As long as, the final temperature of our model 
system stays above T = 0.2, long-time diffusive behaviour 



















5 


* 

-to 




10 ^ 


10 ° 


10 ‘‘' 


10-2 


10 -° 



FIG. 5. (color online) Non-Gaussian behaviour of self-displacement fluctuations Gs{r,t*) (t* = 50t) is shown for compact 
and non-compact aggregates in (a) and (b) respectively. Please note the change in microscopic dynamics accompanying the 
structural transition. While dominant Gaussian self-fluctuations end with an exponential tail for compact clusters, displacement 
of a particle from its arbitrary origin is distributed exponentially over space for non-compact morphologies, (c) Temporal 
evolution of dynamic Lindemann parameter 'y(t) measuring the relative displacement fluctuation of particles with respect to 
its neighbours captures another aspect of this dynamical transition. Very slow initial increase in 7 (t) for compact aggregates 
clearly points to the caging of particles by their near-neighbours. On the other hand, initial steep increase followed by significant 
slowing down of 7 (t) at long-time can be attributed to bonding in non-compact structures. 


is observed (FigQa)). Sub-diffusion at long-time appears 
for T < 0.2. We show the MSD for both compact (^ = 0.5) 
and non-compact = 0.8) structures at the lowest tem¬ 
perature T = 0.05 reached by cooling in our simulations. 
Though sub-diffusion is evident for both types of struc¬ 
tures, the exponent n varies. Interestingly, n for non¬ 
compact string-like structures is the same as the expo¬ 
nent for single-file diffusion [46]. One important feature 
revealed from the MSD data is the prominent intermedi¬ 
ate slow relaxation for compact morphologies. This fea¬ 
ture is also observed for non-compact ones, though with 
less prominence. This intermediate relaxation is indica¬ 
tive of heterogeneous dynamics which we shall probe in 
the next section. 

Here, we continue investigating average dynamics fur¬ 
ther to study the overlap between two arbitrarily chosen 
instantaneous configurations to probe the degree of struc¬ 
tural relaxation within the system. Quantification of the 
overlap gaiiH] is computed as follows: 

N 

rj{t) = w{\Yi{t) - rj(O)l) ( 10 ) 

hi=i 

where i, j are particle indices and w{r) is a window func¬ 
tion which is unity if |ri—r 2 | < a and zero otherwise. 
Allowing a particle to relax within certain threshold dis¬ 
tance a from its original position, this quantity tells us 
whether a particle is knocked out by its own copy or any 
other particle at later time. We have chosen a = 0.03cr as 
this is the average distance particles traverse in unit time 
r. From the MSD data, we can easily identify this length 
scale and the corresponding time scale fall within the in¬ 
termediate slow relaxation regime for all low temperature 


morphologies. In FigQb), we show that configurations 
at T = 0.05 become uncorrelated in a scale-free manner 
from their arbitrarily chosen initial conformations after 
certain residence time r^. We will return to the depen¬ 
dence of Tr with different aggregate shapes later. Here, 
we note that this power-law behaviour for structural re¬ 
laxation changes drastically when the final temperature 
of the system goes above T = 0.2. For these high tem¬ 
perature systems, the threshold a = 0.3cr is again the 
average extent of particle displacement in unit time r 
and the system goes from ballistic to diffusive after this 
time. Considering the temperature dependence of self¬ 
overlap of configurations and MSD of constituent parti¬ 
cles together, we infer that the underlying phase separa¬ 
tion process affects the average dynamics of the system as 
the system crosses gas-liquid critical temperature during 
cooling. However, local dynamics is governed by thermal 
fluctuations set by the final constant temperature. 


C. Local dynamics: self and relative displacement 
fluctuations 

The statistics of individual displacement fluctuations 
for different local geometries is presented in this section. 
The probability of finding a normally diffusing particle 
is given by a Gaussian distribution with variance related 
to its long-time diffusivity. However, this distribution is 
often complemented by a fat exponential tail for systems 
exhibiting anomalous diffusion such as ours. We deter¬ 
mine the probability of finding a particle at a distance 
r(t) from its arbitrarily chosen previous position r(0) by 
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computing the self-part of van Hove function [49] , 

Gs{r,t) = 1/N'^{5{r - |ri(t) - rj(0)|)) (11) 

i=l 

Interestingly, a broad range of microscopic dynamics are 
captured by Gs (r, t) for aggregates across their morpho¬ 
logical hierarchy. At the end of intermediate slow relax¬ 
ation (t* =50r), Gs{r^V^) is dominantly Gaussian with 
small exponential tail (Figj^a)) for compact clusters. On 
the contrary, Gaussian part of Gs(r, t) is recessive and 
buried within the long exponential part (FigJ^b)) for 
non-compact aggregates. These two distinct parts im¬ 
mediately point to at least two different corresponding 
diffusivities m other than the average diffusivity and 
thus, the heterogeneous dynamics in the system. Most 
importantly, by tuning the competing interactions, we 
are able to capture the whole range of spatio-temporal 
heterogeneity for aggregates. 

Next, we identify different physical mechanisms re¬ 
sponsible for the limiting cases of such heterogeneity us¬ 
ing the evolution of relative displacement fluctuations. 
This is captured by calculating the dynamic Lindemann 
parameter EB, 

lit) = ((Aui(t) - Auj{t)f)/2a^ (12) 

as the difference between displacement Au(t) = r(t)—r(0) 
of Tth particle and its j-th neighbours. The particles in 
compact aggregates (^ = 0.5, 0.6) show negligible or very 
little displacement with respect to their neighbours un¬ 
til time ^ r (Figj^c)). This leads us to conclude that 
the particles, in this case, are caged by their neighbours. 
After certain time, the particles get out of the cage and 
can move freely till they get into another cage and 7(t) 
increases monotonically with time. However, from visual 
inspection, we expect this to happen more frequently for 
the particles at the edge of a cluster and less frequently 
for those at the centre of compact crystalline aggregates. 
A reversal of this feature is observed for non-compact 
aggregates. Large repulsion = 0.7,0.8) responsible 
for these chain-like structures hinders the particles from 
coming close to each other and 7(t) shows a sharp initial 
increase. Once they are able to decrease their relative 
distance due to thermal fluctuations, strong attraction 
binds them and shows considerable decrease in the 
slope. We name this feature as bonding. In essence, as 
the aggregate morphology changes from compact to non- 
compact^ the dynamics of constituent particles show tran¬ 
sition from caging to bonding. 


V. DISCUSSION 

The mobility of particles forming aggregates becomes 
extremely sluggish due to the nominal thermal fluctua¬ 
tions present in the system at low (T < 0.2) temperatures. 
The microscopic dynamics is heterogeneous following the 
restricted and random spatial extent of aggregates. Fol¬ 


lowing the individual particle trajectories in different ag¬ 
gregate morphologies and the evolution of their fluctu¬ 
ations, we have identified the fundamental mechanisms 
for this continuous transition in heterogeneous dynam¬ 
ics. Under strong influence of short-range attraction, 
the movement of a particle is restricted within a cage 
bounded by its neighbours and particles align themselves 
in local hexagonal arrangement, energetically favoured in 
2D. Consequently, the probability of finding a particle at 
any instant will be a Gaussian with variance equal to the 
typical cage size determined by the time window between 
the chosen instant and an arbitrary time origin. The 
probability will decay exponentially outside this typical 
size due to the overall rearrangements happening within 
the system at that instant. With increasing repulsion, 
the degree of geometric frustration due to competition 
between interactions increases and the hexagonal sym¬ 
metry becomes unfavourable. An emergent anisotropy 
originating from the complex interplay between enhanced 
frustration and thermal fluctuations allows only sidewise 
bonding of particles. As a result, particles assume finite- 
size string-like non-compact structures and retain their 
shape as the probability of finding the constituent par¬ 
ticles becomes exponentially smaller with increasing dis¬ 
tance from their original position. 

Following this physical picture, it is only natural to ex¬ 
pect that particles within a crystalline compact arrange¬ 
ment stay longer at or in close proximity to their original 
position compared to those in non-compact conforma¬ 
tions. Our computation of the self-overlap function, in 
fact, supports this expectation (FigQb)). The long-time 
algebraic tail of self-overlap can be explained in terms of 
a survival function of the form {rr/t)^ where is the 
minimum residence time of a particle and 7 is the tail 
or Pareto index [52]. The typical residence time is, in¬ 
deed, large for compact clusters and drops down at least 
by amount r for non-compact conformations. Here, we 
mention that these typical time scales are specific to the 
local geometry and not to the specific values of param¬ 
eter f. We recall that aggregates of statistically similar 
morphology can be obtained for various combinations of 
attraction range (fixed by a), repulsion strength (A) and 
screening length {f). The morphological hierarchy can 
be classified by using a dimensionless parameter A [39] 
which encodes all three control parameters {a. A, of 
the effective interaction. As A faithfully describes the lo¬ 
cal geometry of aggregates, we plot corresponding to 
these geometries with respect to A in FigJ^a). A clear 
transition in is observed at A^2. which coincides with 
non-compact (A>2.0) to compact (1.5<A<2.0) transi¬ 
tion. However, A itself is a function of temperature. For 
higher values of final temperature (which means the cool¬ 
ing was stopped at higher temperature), thermal fluctua¬ 
tions compete with geometric frustration and finally take 
over at a temperature above T = 0.2 and the system be¬ 
haves like an equilibrium liquid. 

Finally, we proceed further to relate the average dif¬ 
fusivity with the structural randomness of the system. 
Gonsidering the diffusion in equilibrium simple liquids is 
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FIG. 6. (color online) (a) Temperature dependence of re¬ 
duced diffusivity D* is plotted against the dimensionless de¬ 
scriptor A of aggregate morphologies. Please note that A is 
also temperature dependent, it assumes lower values with in¬ 
creasing T. For low temperatures, T < 0.2, D*-K relation 
always shows an inflection point around the structural tran¬ 
sition. At and beyond T = 0.2, D* becomes constant for 
all A as the system behaves like an equilibrium liquid, (b) 
These observations can be understood more clearly in terms 
of residence time of particles Ta as a function of A. For ex¬ 
ample, Ta marks a clear dynamical transition with respect 
to A at T = 0.05. This transition goes away with increasing 
temperature, (c) D* can further be related to configurational 
randomness of the system measured by two-body excess en¬ 
tropy 82 - In fact, D*-S 2 follows the universal relation {see 
text) for small values of aS '2 or high randomness. The compact 
aggregates at low temperatures having local crystalline order 
naturally deviate from this universality. 

governed by the Enskog process [53], it is possible to re¬ 
late the diffusivity of the system to the configurational 
phase space volume available to the same, that is, 

D = Ca^Te^^ (13) 

where F is the typical collision frequency and £' 2 , the 
two-body excess entropy [54| , is an ensemble independent 
measure of structural randomness of the system [55]. C 
is a system dependent dimensionless non-zero constant. 
Such relation was phenomenologically proposed and es¬ 
tablished for simple equilibrium model systems [56] [57] . 
Later, it was experimentally verified for colloidal systems 
[58] and a wide range of liquid metals [59] . For our model 
aggregates, non-compact clusters posses higher degree of 


randomness {S 2 ^ 0) closer to liquids. S 2 assumes large 
negative values as soon as local crystalline ordering ap¬ 
pears within the system as in compact aggregates [39] . 
Since for all these morphologies, the microscopic dynam¬ 
ics of our model system at low temperature is non-ergodic 
and much slower than simple liquids, an Enskog process 
cannot be expected. We define a dimensionless diffusiv¬ 
ity D* = D^/Dq where = limt-^(r^(t)) is the so- 
called long time diffusivity and Do = limt-^{r‘^{t)) repre¬ 
sents particle motion in short-time ballistic motion when 
the particles have yet to experience the presence of their 
neighbours. Plotting 1)* against S 2 (Eig j^b )), an expo¬ 
nential relation, D* ^ (similar to eqfl^ is observed 
for different temperatures and small values of 82 - We 
note that small S 2 {< 1 ) refers to equilibrium liquids for 
T > 0.2 and non-equilibrium non-compact structures for 
T <0.2. This observation is surprising and implies that 
the Enskog process is not limited to equilibrium but re¬ 
lated to the morphological randomness. As the system 
adapts specific crystalline symmetry this feature breaks 
down and I)* becomes nearly independent of 82 for the 
compact aggregates. With the possibility of a fundamen¬ 
tal description of microscopic dynamics in terms of the 
degree of configurational randomness, this observation 
surely requires further study. 

To conclude, aggregates are suitable candidates to re¬ 
late non-ergodic dynamics to the corresponding local 
morphology. Geometric frustration inherent to the com¬ 
peting short-range attraction and long-range repulsion 
leads to non-ergodic stationary aggregates. The average 
dynamics bearing the effects of phase separation at suf¬ 
ficiently low temperature shows long time sub-diffusive 
behaviour following intermediate slow relaxation. Both 
structural randomness of aggregates and the heteroge¬ 
neous micro-dynamics of constituent particles can be 
changed seamlessly by tuning the interactions for fixed 
thermodynamic conditions. Specifically, the particles in 
compact aggregates experience caging from their nearest 
neighbours and individual dynamics are close to that of 
particles in glassy systems. On the other hand, bonding 
of particles is prevalent in non-compact clusters and the 
individual dynamics resembles that of particles moving 
in a network. Eurther correspondence between dynamics 
and morphology of aggregates has been established by 
relating diffusivity and two-body excess entropy. In fu¬ 
ture, we plan to investigate the effect of cooling rate, age¬ 
ing time and external mechanical fields on these pattern 
forming systems to understand the role of local geometry 
in systems with unusual dynamical responses. 
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